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Abstract 



Collective P- and T- odd moments produced by parity and time invariance 

> . 

violating forces in reflection asymmetric nuclei are considered. The enhanced 

o 

(N 



collective Schiff, electric dipole and octupole moments appear due to the mix- 



\q \ ing of rotational levels of opposite parity. These moments can exceed single- 

particle moments by more than two orders of magnitude. The enhancement 
is due to the collective nature of the intrinsic moments and the small energy 
separation between members of parity doublets. In turn these nuclear mo- 



ments induce enhanced T- and P- odd effects in atoms and molecules. First a 



simple estimate is given and then a detailed theoretical treatment of the col- 

lective T-, P- odd electric moments in reflection asymmetric, odd-mass nuclei 

is presented and various corrections evaluated. Calculations are performed 

for octupole deformed long-lived odd-mass isotopes of Rn, Fr, Ra, Ac and 

Pa and the corresponding atoms. Experiments with such atoms may improve 

substantially the limits on time reversal violation. 
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I. INTRODUCTION 



In 1964 Christenson et al. have discovered that CP is violated in the decay of neutral 
K-mesons. As one expects the CPT theorem to be valid this discovery implied that time 
reversal (T) is violated in the observed decays of the Kaon. This fact has immediately lead 
to the search of CP or T- violation in other systems. Since the 60' s many attempts have been 
made to observe T-violation in systems different from the Kaons. Time reversal violation 
has not been observed so far but upper limits for T-conservation have been established. The 
search of time reversal violation encompasses a large variety of physical systems and involves 
many methods. One of the more widely used methods involves the search for static T-, P- 
odd electromagnetic moments, it is moments that would be absent if the Hamiltonian of 
the system is even under time reversal and reflection. Such moments include the electric 
dipole moment (EDM), the electric octupole, the magnetic quadrupole etc. Early on, with 
the discovery of CP violation, attempts were made to measure the electric dipole moment 
of the neutron and at present significant upper limits exist on the existence of such moment 
@||. The neutron was not the only system in which attempts were made to find a static 
electric dipole moment. Experiments with atoms and molecules were performed in which 
upper limits for electric dipole moments of the respective systems were established. In fact 
the recent measurements of dipole moments of Hg and Xe atoms M and TIF molecule || 
have established upper limits for time reversal violating nucleon-nucleon and quark-quark 
interactions that are of the same order (or maybe even exceed) the limits obtained in the 
measurement of the neutron dipole moment. 

The existence of a static atomic dipole moment may be due to the following three reasons: 
a) the possible existence of a dipole moment of the electron (the best limits on the electron 
EDM were obtained in Tl atom EDM measurements in Ref. 0, see also Cs measurements 
in Ref. J/J) b) time reversal violation in the electron-nucleon interaction, thus in the lepton- 
hadron interactions, c ) the possible existence of a static T-odd, P-odd moment of the nucleus 
arising from the time reversal violating component of the hadron-hadron interactions. The 
recent experiments with Hg gave the best limit on this interaction. This possibility will also 
be the subject of this work. 

In a recent paper M we put forward a suggestion that rotating nuclei that have static 
octupole deformations when viewed in their intrinsic (body) frame of reference will have 
enhanced T-odd, P-odd (for short T-, P- odd) moments if a time reversal and parity violating 
interaction is present in the nuclear Hamiltonian. In the intrinsic frame the nucleus with 
an octupole deformation has large octupole, Schiff and dipole moments. An orientation 
of these moments is connected to a nuclear axis n (e.g. the dipole moment is d = dn). 
In a stationary rotational state the mean orientation of the axis vanishes ((n) = 0) since 
the only possible correlation (n) oc / violates time reversal invariance and parity (here I is 
the total angular momentum of the system). Therefore, the mean values of electric dipole, 
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octupole and Schiff moments vanishes in laboratory frame if there is no T-, P- violation. In 
the nuclei with the octupole deformation and non-zero intrinsic angular momentum there 
are doublets of rotational states of opposite parity with the same angular momentum I (in 
molecular physics this phenomenon is called A - doubling). A T-, P- odd interaction mixes 
these rotational levels. As a result the nuclear axis becomes oriented along the total angular 
momentum, (n) oc al where a is the mixing coefficient. Due to this orientation of the 
nuclear axis by a T-, P- odd interaction the mean values of the T-, P- odd moments are not 
zero in the laboratory frame, e.g. (d) =d(n) oc adl. 

We find two basic enhancement factors in this mechanism: firstly, in the intrinsic frame 
the nucleus with an octupole deformation will have large octupole, Schiff, or dipole electric 
moments because a large number of nucleons will contribute to the moments, and secondly 
due to the appearance of closely spaced parity doublets in the spectrum of the nucleus with 
octupole deformation. It is not only that the spacing between the members of the doublets 
is small but also (T-, P- odd) interaction will mix well such two states. The enhanced 
nuclear Schiff moments that result in such nucleus with a reflection asymmetric shape will 
induce ~ 1000 times enhanced atomic electric dipoles, and measurements performed with 
such atoms may improve upper limits for time reversal violation. 

It is the aim of the present work to examine in detail the consequences of the intrinsic 
reflection asymmetry on the T-, P- odd electromagnetic moments in nuclei produced by T-, 
P- odd components in the nuclear force and on the induced T-, P- odd moments in the 
corresponding atoms. In this paper we extend the work in Ref. || attempting to provide an 
improved and more detailed theory of the nuclei with asymmetric shapes and of the resulting 
T-, P- odd moments if parity time reversal is violated to some degree in the nuclear force. 

Present experimental studies of nuclei in the actinide region (Z around 88 and N around 
134) indicate that these nuclei posses octupole shapes in the ground state (g.s.) PHT^B- 
In these nuclei near the Fermi energy, orbital pairs are coupled strongly by the octupole- 
octupole part of the effective nuclear interaction. The existence of octupole deformations in 
the actinide nuclei is manifested in the existence of parity doublet states and parity doublet 
bands. The El and E3 transitions between these states are relatively strong, of the order of 
a Weisskopf unit. These experimental findings are supported by theoretical studies. Some 
isotopes of Rn, Fr, Ra, Ac, Th, and Pa in the 218<A<230 region are predicted theoretically 
to be reflection asymmetric in the g.s. The phenomena of octupole instability are observed 
and described theoretically also around Z ~ 56, N ~ 88. Several isotopes of Ba, Ce, 
Nd, and Sm in the A=140-152 region are known to be octupole-soft and develop reflection 
asymmetric shapes at higher spins but no experimental data at present can confirm such 
shapes in the g.s. |TT|JT2| . 

We will present in this work results for relatively long-lived neutron-odd isotopes of Rn, 
Ra, and proton-odd isotopes of Fr, Ac, and Pa for which there is theoretical and in most 
cases experimental evidence of reflection-asymmetric intrinsic shapes in the g.s. 
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We use the results of the Nilsson- Strutinsky mean field calculations Jl3| and employ the 
Leander et al. particle plus core model to calculate the wave functions, energy splittings 
and T-, P- odd interaction matrix elements between members of the parity doublets. In the 
calculation of the intrinsic dipole, Schiff and octupole electric moments we use a two-liquid 
drop model. Strutinsky corrections and corrections due to pairing are taken into account in 
this work. 

In the past one conjectured that T-, P- odd moments will be enhanced in the cases 
when close to the g.s. there is, for whatever reason, a level with the same spin and opposite 



parity p^-l8[. The connection between enhanced El transitions between these levels and 
the nuclear electric dipole moment was made in the work of Haxton and Henley [jT7|| . At 
that time, nuclei which possess an intrinsic reflection asymmetry and quadrupole deformed 
nuclei without such an asymmetry but in which an accidentally close "parity doublet" exists, 
were treated on an equal footing ||17|| . The fact that in some nuclei systematically enhanced 
El matrix elements are the direct consequence of intrinsic reflection asymmetry was realized 
and investigated extensively by Leander et al. several years later ||19| . 

In the work []TH] the Schiff moment induced by nuclear T-, P- odd forces was introduced 
in the presently used form and calculations for Xe, Hg, Tl and other interesting cases were 
done. The calculations of the atomic electric dipole moments induced by the nuclear T-, P- 
odd Schiff moments were also presented (we should note that similar considerations for T-, 
P- odd effects in molecules and atoms induced by the proton electric dipole moment were 
applied in J2D| and Q). 

In our earlier paper || for the first time the connection was made between the collective 
T-, P- odd electric moments in the intrinsic frame of reference in reflection-asymmetric 
nuclei and these moments in the laboratory frame. It is in this case, due to the collective 
nature of the intrinsic moments and the nearly identical intrinsic structure of the parity 
doublets one that hopes the electric T-, P- odd moments will be maximal. 

After the present introduction in Sec. fTI] we define the T-, P- odd moments, including 
the Schiff moment. In Sec. |TT| we present a simple expression for a T-, P- odd moment 
in the case of a deformed rotating nucleus in the presence of a T-, P- odd interaction. In 
the same section we bring a simple schematic estimate of a Schiff moment in an octupole 
deformed nucleus. The first part of Sec. [TV] deals with the calculation of the dipole, Schiff 
and octupole intrinsic moments in a two-fluid liquid drop model. In the second part of 
this section we present the particle+core model and describe the calculation of the T-, P- 
odd matrix elements and mixing amplitudes. In Sec. the numerical results are presented 
for each of the nuclei and at the end of this section results are given for the atomic dipole 
moments. In the last section (Sec. |VI|) a summary is presented. 
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II. ATOMIC ELECTRIC DIPOLE AND NUCLEAR T-, P- ODD MOMENTS 



We start from the electrostatic potential of a nucleus screened by the electrons of the 
atom. If the nucleus has a T-, P- odd dipole moment the dipole term in the potential 
vanishes in accordance with the Purcell-Ramsey-Schiff theorem [0,^] (see the derivation 
in Appendix |V|): 

<p(r) = I #V»- + i(rfv) / • (i) 



\R-r\ Z K 'J \R 

Here p(r) is the nuclear charge density, J p(r)d 3 r = Z, and d = J erp(r)d 3 r is the electric 
dipole moment (EDM) of the nucleus. The first term in this expression is usual electrostatic 
nuclear potential, and the second term is a result of the electron screening effect. The 
multipole expansion of <p(R) contains both T-, P- even and T-, P- odd terms. We consider 
here only the latter. The dipole part in Eq. ([!]) is canceled by the second term in this 
equation: 

- J e(rvi)p(r)d 3 r + ^(dV)± J p(r)d 3 r = . (2) 

The next term is the electric quadrupole which is T-, P- even, thus the first non zero T-, P- 
odd term is: 



If 1 

^ (3) = ~7 / ep(r)r a r /3 r 7 rf 3 rV Q V /3 V 7 - 
J K 

+ 2^ (dV)VaV/3 | / P^ r >^ r ■ ( 3 ) 

Here r a rpr^ is a reducible tensor. After separation of the trace there will be terms which 
will contain a vector S and a rank 3 tensor Q a /3-y (see e.g. |p~8l ): 

^(3) _ (3) , (3) 

V — ^Schiff ^ 'foctupole ' 

As = " 5VA ^ = 4nSV5(R) , 

^ictupolc = -g<5«/3 7 V « V /3V 7 — , (4) 



where 



S = — 
10 

is the Schiff moment (SM) and 

1 



— (J ep(r)r 2 rd 3 r - -d— J p{r)r 2 d z r^j (5) 



Q a /3j =J ep(r)[r a r0r 7 - |(5 a/3 r 7 + 5 Pl r a + 8 ai r )]d 3 r 
Q zzz = Iq 3 = IJyJ ep(r)r 3 Y 30 d 3 r (6) 
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is the electric octupole moment. Because the intrinsic dipole moment of the nucleus appears 
in second order in the nuclear deformation (see Eq. fl34|)) a correction to the octupole field 
which arises from the non-spherical part of the density in the screening term in Eqs. (0,0) 
can be neglected. Indeed, this correction is at least third order in the nuclear deformation. 

In the absence of T- and P- violating interactions the electric dipole moment of an atom 
is equal to zero. The interaction between atomic electrons and the T-, P- odd part of the 
electrostatic nuclear potential in Eq. (f4|) will mix atomic states of the opposite parity and 
thus generate an atomic electric dipole moment: 

D z = -e(i/>\r s \i/>) = -2e}_^ , (7) 

\k 2 ) _ ^2 

where if) denotes the perturbed atomic wave function, \ki) = \ki,ji,ju) is the unperturbed 
electron ground state and {^2)} is the set of opposite parity states with which \ki) is mixed 
due the perturbation —e(p^. 

The most accurate measurements of atomic and molecular T-, P- odd electric dipole 
moments have been done in the atoms Xe and Hg with zero electron angular momentum, 
ji = 0. Examining Eq. (|7|) it is easy to demonstrate that in such atoms nuclear electric 
octupole (as well as another T-, P- odd moment, magnetic quadrupole) cannot generate an 
atomic electric dipole. Indeed, according to the triangle rule for the addition of angular 
momenta, (&i|r 2 |/c 2 ) can only have a nonzero value if \ji — j 2 \ < 1 < ji + j 2 - Similarly, for 
(ki I <^octupoie I ^2) to be nonzero, we must have \ji — j'2 1 < 3 < ji + j 2 . This implies that the 
following conditions need to be satisfied for the dipole moment to be nonzero: 

|ji - J2I < 1 and ji + j 2 > 3 . (8) 

The lowest pair of values that satisfies this condition is ji = 3/2 and j 2 — 3/2 for the states 
with one or odd number of electrons outside the closed subshells, and j\ — 1, j 2 = 2 for states 
with even number of electrons. In the case of the magnetic quadrupole one needs j\ + j'2 > 2, 
i.e. the lowest pair is ji = 1/2, j 2 = 3/2. Hence, the nuclear electric octupole and magnetic 
quadrupole moments cannot contribute to the atomic electric dipole moment of ji = states. 
Therefore, we mostly center our considerations on the Schiff moment. Note also, that the 
T-, P- odd part of the electrostatic nuclear potential in Eq. (|) is concentrated mainly inside 
the nucleus. As a result the induced atomic electric dipole moment is proportional to the 
density of external electrons at the nucleus (more accurately, to the gradient of this density) 
which rapidly increases with the nuclear charge Z. This is why in general heavy atoms and 
nuclei are favored in the studies of T-, P- odd moments. 



III. SIMPLE ANALYTICAL ESTIMATE OF THE T-, P-ODD MOMENTS 
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A. T-, P- odd moments and rotational doublets 



If a deformed nucleus in the intrinsic (body-fixed) frame of reference is reflection asym- 
metric, it can have collective T-, P- odd moments. As a consequence of this reflection 
asymmetry, rotational doublets appear in the laboratory system. Without T-, P- violating 
forces a T-, P- odd moment vanishes exactly in the laboratory. A T-, P- odd interaction 
however may reveal such intrinsic T-, P- odd moments in the laboratory frame. Consider 
a nearly degenerate rotational parity doublet in the case of an axially symmetric nucleus. 
The wave functions of the members of the doublet are written as 



^ ± = -=(\IMK) ± \IM- K)) . (9) 
v2 

Here I is the nuclear spin, M = I z and K = In, where n is a unit vector along the nuclear 
axis. 

The intrinsic dipole and Schiff moments are directed along n: 

S'intr — •S'intr/l . (10) 

For these good parity states {^ ± \In\^/ ± ) = because K and —K have equal probabilities 
and this means that there is no average orientation of the nuclear axis in the laboratory frame 
((^ r± |n|^ r± )= 0). This is a consequence of time invariance and parity conservation since the 
correlation In is T-, P-odd. As a result of (\E' ± |n|\l/ ± )= 0, the mean value of the T-, P- 
odd moments (whose orientation is determined by the direction of the nuclear axis) is zero 
in the laboratory frame. 

AT-, P- odd interaction Htp will mix the members of the doublet. The admixed wave 
function of the predominantly positive parity member of the doublet will be ^ = + a^>~ 
or 

¥ = -^=((1 + a)\IMK) + (1 - a)\IM - K)) , (11) 



where a is the T-, P- odd admixture 

and E + — E~ is the energy splitting between the members of the parity doublet. A similar 
expression is obtained for the negative parity member of the doublet. 
In the T-, P- admixed state 

(■9\In\9) = (^\K\^) = 2aK , (13) 

i.e. the nuclear axis n is oriented along the nuclear spin J: 
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mnM = 2a irrryi ■ (14) 

Therefore in the laboratory system the electric dipole and Schiff moments obtain nonzero 
average values. For example in the ground state (g.s.) usually M = K = I and 

(*\S,\V) = 2aj^S intI . (15) 



B. Estimate of the Schiff moment 

Let us present a simple estimate of the nuclear Schiff moment in the case there is a T-, 
P- odd interaction and the nucleus possess a quadrupole and octupole deformation. One 
needs to calculate the intrinsic moment and the mixing coefficient a. We consider an axially 
symmetric nucleus which has a sharp surface and constant nucleon density p t . The surface 
of the nucleus has the form 

R = R (l + &Y 10 + (3 2 Y 20 + (3 3 Y 30 ) . (16) 



The f3\ deformation is introduced p5,26,24| to retain the center of mass at z = 0: 



/ 



Thus 



4-TT 9 

zp t d 3 r ~ Bo(y T A + -^mlhth) • (17) 



ft= -VsiH ftft ' (18) 



The proton density can be approximated by constant p = 3Z/AttRq. The intrinsic electric 
dipole moment in this approximation is zero, and only the first term in Eq. (|5p contributes 
to the intrinsic SM. Using Eqs. (|,|1^,|1^) one gets: 

S intr = ^e[ pr 2 zd 3 r = eZI^—^fcfc . (19) 
1U J 2()7rv35 



207TV35 

Using then the deformation parameters p 2 — 0.12 and (3$ = 0.1, Rq = 1.2 x A 1 / 3 fm, 
A = 230, and Z = 88 we obtain : S intr = 10.4e fm 3 . 

The rotational states of an odd-mass nucleus can be written in terms of intrinsic states 

as: 

\IM ±K)= (^±1) l ' 2 D{ I±K ^ 9, ^)^(r')x (A) , (20) 

where D^^^, 9, ip) is a Wigner D-function, is the wave function of the quadrupole and 
octupole deformed (reflection asymmetric) nuclear core in the intrinsic frame, and y?±^-(r') is 
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the wave function of the unpaired nucleon in the intrinsic frame, with an angular momentum 
projection of ±K on the z' axis. 

Let us now present an order of magnitude estimate of the mixing coefficient a which is 
needed to find the magnitude of the Schiff moment in the laboratory frame (see Eq. (|T5D). 
K = In and Htp are both T-, P-odd pseudoscalars. Therefore, (p^\Htp\p^) oc K 
and so (^-k\Htp\v^k) = — (^+k\Htp\^+k) (this fact can be easily supported by model 
calculations). Using this fact and Eqs. (|20| ) and (|) we get (^/~\H TP \^ + } = (p^\H TP \p^). 
If single-particle wave function p^ was a good parity state this matrix element would be 
zero. However, due to the perturbation caused by the static octupole deformation of the 
nucleus (V3), it is a combination of the opposite parity spherical orbitals 4>i,+k and 4>2,+k 
(for example, S1/2 and pi/2, or p 3 / 2 and (I3/2, etc.): 

tfi+K = <Pi,+k + I<p2,+K , (21) 
7= ( ^y , (22) 

— -&2 

4> lt+K = Ki(r')VL j)l)+K (6' ,p') j 
= Mr')Sl jJj+K (9',<p') 

= -n 2 (r')(<Tr')n^ +K (V',<P') , (23) 

where I = 2j — I. (Of course there will also be an admixture of other opposite parity states 
having different values of j . We neglect these states for simplicity. ) The mixing coefficient 
7 is proportional to the parameter of the octupole deformation /5 3 . Using Eq. (|22D one can 
check that numerically 7 ~ /?3. Therefore, we have 

_ W^HtpW+k) (<Pi,+k\H T p\<P2,+k) (oa s 

a - E+-E- ~ m E^W- • (24) 

Finally, we 
must estimate the matrix element between spherical orbitals {4>i,+k\Htp\4>2,+k) ■ The T- 
and P- odd interaction between a nonrelativistic unpaired nucleon and the nuclear core can 



be described by the following effective Hamiltonian (see e.g. Refs. p8|j27|j28[] ) : 



G 

Htp = 27— =-crVp t , (25) 
2\/2m 

where er is twice the spin operator for this nucleon, p t is the nucleon density of the nuclear 
core, G — 1.0 x 10 _5 /m 2 is the Fermi constant, m is the mass of the nucleon and 77 is a 
dimensionless constant that describes the strength of the interaction. 

Using Eq. fl2"o| ) for the form of Htp and ptir 1 ) = 6{r' — _R )/(f 7r r o 3 ) (where 9 is the step 
function) we get 
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H TP = -V 8n j£ 3 (<r*')W - Ro) ■ (26) 
Here ro = 1.2 fm is the internucleon distance. Using Eq. (^) and (err) 2 = 1 gives 
\4>\,+k\Htp\4>2,+k) 

= ^fe K > (WW ^ leV ' (27) 

where we have used K 1 {R Q )n 2 {Ro) ~ l-4/i? 3 H- Using \E + - E~\ = 50 keV, 
(3 3 = 0.1 (see data in the Tables), and Eqs. ([24]) and (p7] ) gives (for A = 230) |a| ~ 
2PsA~ 1 / 3 r] eV/\E + — E~\ ~7x 10~ 7 7/. This provides the following estimate for the collec- 
tive Schiff moment in the laboratory frame: 

S ~ aS intr ~ 0.05e/3 2 /33 2 ZA 2 / 3 r/ (r ) 3 eV/|£+ - 

~ 700 x 10~ 8 e fm 3 . (28) 

We see that the collective Schiff moment is about 500 times larger than the Schiff moment due 
to the unpaired nucleon in spherical nucleus, S ~ 1.5 x 10~ 8 ?7e fm 3 [18|,|3(|. Note however 



the strong dependence of the collective Schiff moment on the deformation parameters. 

We should remark here that this is a schematic calculation and more detailed and realistic 
calculations are presented in sections |V| and |V|. Nevertheless the present estimate represents 



the essence of this theory and the order of magnitude estimate agrees with the detailed 
calculations. 



IV. NUCLEAR MODELS OF THE T-, P-ODD MOMENTS 

A. Nuclear shape and intrinsic moments 

In this paper we consider the moments of heavy deformed nuclei in the ground states. 
The main contribution to the electric moments comes from the even-even core which is well 



described by the two-fluid liquid drop model, see e.g. Refs. pT| , p9| , p4| . The surface of an 
axially symmetric deformed nucleus is 

R = c v (f3)R (l + J2PiY l0 ) , (29) 
1=1 

where c v = 1 -== ensures the volume conservation and Rq = tqA 1 ^ 3 . For the sake 

V47T; =1 

of brevity and because the nuclear deformations we deal with are relatively small fa < 0.2, 
/?3,4 < 0.1 we will put in our discussion the coefficient c v ((3) = 1. (In the actual calculations 
Ro is replaced by c v ((3)Rq.) 
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If the nucleus is reflection-asymmetric then the (5\ deformation parameter is needed to 
keep the center of mass fixed at z — 0, i.e. / zd 3 r = 0. In lowest order in nuclear deformations 

m 



-3#E 

V Air ^ 

v 1=2 



(z + i)AA+i 

(2/ + l)(2/ + 3) 



(30) 



Due to the Coulomb force protons and neutrons are differently distributed over the nuclear 
volume. From the requirement of a minimum in the energy [j3~l),|19|1 



P P (r) ~ Pn{r) 
p P (r) + p n (r) 



(31) 



where p p = p and p n are the proton and neutron densities, Vc ou i( r ) is the Coulomb potential 
created by p p (r) and C is the volume symmetry-energy coefficient of the liquid-drop model. 



To lowest order 19 



Po _ Po e 2 Z 
2 8 CR 

3 1 / r \2 „ 

(id + S 



x 



2 ^i?o i =1 2/ + 



(32) 



where po = 3A/(4tiRq). The coefficient C is not known very accurately, its value for nuclei 
studied here is in the range 20-35 MeV [ 19U26 ,|32]. Note that requiring / p(r)d 3 r = Z one 



has in lowest order [31 



9, V 



3 e 2 Z 

TocK 



) 



(33) 



We compute the intrinsic Schiff moment by substituting the density in Eq. (|3"2]) into Eq. 
(^). Because of the relative shift of protons versus neutrons the nucleus in the intrinsic frame 
has a dipole moment as calculated in the past |^3lR3],O^BJ and given by 



d 



intr 



eAZ 



(Z 2 -l)(8Z + 9) 
C40tt^ [(2/ + l)(2/ + 3)] 3 / 2 



E 



i+i 



(34) 



A more detailed treatment of the intrinsic dipole moment includes also the neutron skin 
effect which reduces di n t r somewhat PBl|32[. To discuss corrections to the intrinsic SM it is 



convenient to decompose it into two terms 



q _ o(i) I c( 2 ) 

'Jintr °intr ' ^intr 



(35) 



The first term includes only contribution from the constant part of the density, p(r = 0) in 
Eq. fl3"2p. In lowest order in deformation it equals to 
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°40vrV RoC8>£z J(2Z+l)(2i + 3) 



As observed above this contribution comes from the first term in Eq. (Q) only. The second 
term is due to the Coulomb redistribution of the proton density and stems from the last two 
terms in the brackets in Eq. ([32]). A simple derivation gives 



cCO _ p4R 3 3 e 2 Z 29 v (/ + l)AA + i 
^ ~ 407f RoC 28Q 2, / (2/ + 1)(2/ + 3) • ^ 



This term gives about 10% contribution for nuclei with Z ~ 90. Using Eq. (|33|) one can 
approximate the intrinsic SM as 

S intI ~ eZR 3 — V ; (/ + . (38) 

°207r^ / (2/ + 1)(2/ + 3) 



The contribution of the j3^j3^ term in both di ntT and Si ntr , for the nuclei studied here, is of 
about the same size as the contribution of the term. 
The expression for the intrinsic octupole moment is ||24|,|33 



Q3 intr = ^^(/?3 + 2 j-m + T^F&A + -..) • (39) 

2v77T 3 V 7T llv 71 " 

Various nuclear surface corrections to the density of nucleons such as the neutron skin 
are not included in the above equations for the SM. For the intrinsic dipole moment these 
corrections were included in Refs. [^,^] using the droplet model. The corrections for di ntr 
are of the same order as the main term in Eq. (R3I) but have the opposite sign. One can 



conclude from Eqs. (|34] - |3"8"D that since such corrections alter only the term S\ n l T they can 
contribute to the Schiff moment at most 10%. 

Values of the Schiff moment obtained using Eq. ( |3"8"D are about 30% less than given by the 
direct calculation using Eqs. (|5|,|3"i2D (with C ~ 27 MeV corresponding to the value obtained 
in the droplet model which includes the effect of a neutron skin PB"| , |3"2"|j ). This is partly due 
to terms of higher orders in the deformation not included in Eqs. ( |3"3"| - |3"8"D . 

Both the intrinsic SM and dipole moments are second order in nuclear deformation and 
may turn out to be sensitive to details of the proton density distribution. Because of that 
it is important to take into account quantum mechanical corrections to the liquid drop 
model. Such correction can be included through the Strutinsky shell correction method 
.|HB|. In this method the level density is decomposed into a smooth averaged density and 



a remaining part, fluctuating with the shell filling. The corrected expectation value of an 
one-body operator, e.g. the intrinsic Schiff moment (for the similar treatment of the intrinsic 
dipole moment see Refs. (T^rj^|) is written as a sum of "macroscopic" and shell correction 
terms |36| 
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(40) 



As the "macroscopic" part one takes the liquid drop moment, in the case of the Schiff 
moment S^ tl is given in lowest order in $ by Eq. (^) or by approximate Eq. ([38)). The 
shell correction term is given by [|36| , [T9| , |32 



(41) 



where vf are the BCS quasiparticle occupation numbers and - the smoothed single- 
particle occupation numbers, for the state i. The latter determine the averaged level density. 
(Detailed expressions can be found in Refs. |3^J3^| .) 

It is convenient to express the corrections relatively to the S^ tl . The second term in 
Eq. (^|) is proportional to the intrinsic dipole moment rfi ntr . The shell correction to <ii ntr was 
studied in detail in [|TJ,^]. The results show that the shell correction to <ii ntr is of the same 
order of magnitude as the "macroscopic" dipole moment given by Eq. (0). Therefore the 
subsequent correction to the second term of the intrinsic SM in Eq. ([5]) does not exceed 
(5-7)% of the value of S^ tI . 

The shell correction to the deviation of proton and neutron centers of mass for octupole 
deformed actinide nuclei was also investigated in the Refs. |19|,|32|. The deviation of the 



center of mass can be represented formally as a change of the (3\ deformation parameter (see 
Eq. (|TT|)). Note that the intrinsic Schiff moment is in the lowest order proportional to the 
fii for protons (Eqs. (|3l]|, j36| - |3~8|) ) . Therefore the resulting shell correction for the first term of 
the intrinsic SM in Eq. ([5]) is analogous to the one obtained in p^9| , |32| for the proton center 
of mass. 



B. Particle-core model for a reflection-asymmetric nucleus and P-, T- odd mixing of 

parity doublets 

Any nucleus, whether even or odd mass that is reflection asymmetric may possess in- 
trinsic T-, P- odd moments. Such moments will exist in the laboratory frame only if there 
is also a T-, P- odd mixing between g.s. parity doublet levels. If the even-even nucleus is 
perfectly axial- and spin symmetric, a pseudoscalar T-, P- odd operator cannot mix doublet 
states which have identical intrinsic structure P8| . Besides, the energy splitting within par- 



ity doublets is systematically much less in odd (or odd-odd) than in neighboring even-even 
nuclei P^| , ^4]J41]|] . So in odd (or odd-odd) nuclei the mixing should be considerably larger. 



For this reason we consider odd nuclei in what follows. 

We use here the particle-core model for a reflection-asymmetric nucleus [0,0. The T-, 
P- odd as well as P- odd T- even mixing was studied in this model recently |4"T|,|3l| . This 
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model involves two K n = + and K n = members of a parity doublet of the even-even 
core with the energy splitting E c between them. The wave functions of these states x* are 



projections of the reflection-asymmetric core state x A [0 

1 



X 



V2 



(42) 



where P is the core parity operator. The single-particle states ip A are solutions of the 
reflection-asymmetric single-particle plus pairing Hamiltonian H s . p + i^ pa ir (different forms 
of potentials such as folded Yukawa, deformed Woods-Saxon or Nilsson are used in H s _ p 
|14] , |34]^2] , |38|| ) . The wave functions in the model are [|i~4] , j38|1 : 



^ MK 



21+1 



16tt 2 



1/2 



(43) 



where Rzfa) denotes rotation through an angle 7r about the intrinsic 2 axis. The $ p = $ ± 
are particle-core intrinsic states of good parity p. Denoting the good parity particle states 
<j) n we write (in the matrix notation a<p = J2k a fc0fc) 



a+X + <f> + + Kx <t> 
a-X~<P + + b-X + <t> 

The matrix elements of V PT are given by 

(Hi K \v PT \Hj K ) 

= a + b 



(44) 



y PT |<fe> + a_M0^ PT| 



>K 



(45) 



(As already mentioned the pseudoscalar operator V cannot connect states of an even-even 



spin symmetric and axially symmetric core x w |p8 



The expectation value of a T-, P- odd operator O in a T-, P- admixed state <f>f is 
The matrix elements between core states are 

(x + \o\ x -) = (x a \o\ Xa ) ■ 



(46) 



(47) 



One can write the one-body operator O as the sum of core and particle terms O = O core + O p 
and obtain 



(®i\d\®j > = (x A \Ocor C \x A ){a+ ia- j + K i6_ j) 
+ (^ + |Op|07)a+ 4 6_ j + (0t|6 p |0-) a _ ,•&+ 4 . 



(48) 



For z = j this is just the intrinsic moment in the reflection-asymmetric core-particle state 



X A V A = X A ( a( P + + )• For closely spaced doublets a A 



&_ i and 
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(a + ja_ j + b + ib- j) ~ 5y. With respect to the single-particle contribution to a T-, P- odd 
moment in Eqs. (f46| , f48|) note that the admixture of the doublet level is considerably larger 
than the admixtures of the other levels of opposite parity i.e. an ^> Qjy,i 7^ j. Neglecting 
these "off-diagonal" a^ contributions one gets 

(*t\6\*t) « 2a ii (( X j6 CO reix A ) + (v A |Opb A >) , (49a) 

and 

(^ K \o^ \¥+ K ) 

= (IMlO\IM)(IKlO\IK}($t\6 l *=o\®t) , (49b) 

where / is the rank of the operator. For a vector operator such as S (or d) one obtains in 
accordance with Eqs. (|ID||I1| , |I5] ) 



/(j+i; 



2aii((X A |4ore|X A ) + (^JSpW) • ( 50 ) 



The intrinsic SM of the even Z core (x A \S C ore\x A ) * s gi ven by Eq. (flOf). 



We stress again that the essential difference between the reflection-asymmetric (octupole 
deformed) nucleus and the reflection symmetric (quadrupole deformed) one is that the for- 
mer has intrinsic T-, P- odd moments which are essentially collective since they involve 
contributions of the core nucleons. Note that in a nucleus which has a quadrupole deforma- 
tion but no octupole deformation the ground state and its parity "partner" are built from 
a single good parity core state, e.g. \ + 1 an d therefore the coefficients in Eq. (H4T) are then 



a + = 1, b + = 0, a_ = 0, 6_ = 1. So apart from a small core polarization contribution ||T8|J30 
there is no contribution from the core. 

The parity mixing of a reflection asymmetric single particle state (p A = 
J2k a k<fii + E m fr-m0m * s conveniently expressed via the expectation value of the single-particle 
parity operator 7r p [ |l4|j : 

% = (<Pa\Kp\<Pa) = H a l -H h2 m ■ ( 51 ) 

k m 

Using this quantity one can write the admixture coefficient 





<<^ PT IC> 







(1 - 4i) ■ (52) 

Here E\ and E i are the energies of the particle-core states $>f, $7- The an obviously is 
maximal for a doublet built on the strongly parity admixed (|7r p | -C 1) intrinsic state. Note, 
that although the form of Eq. ( p^) is analogous (for 7r pi = 0) to the case when there is a close 
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single-particle level of opposite parity in reflection-symmetric deformed nuclei ||18|| , there is 
an important difference. Namely [|TJ| , if one neglects the spin-orbit interaction and assumes 
that the nuclear density pt{r) and the single-particle potential U(r) have the same form 
then 

F PT ~<rW(r) =i[crp,H^} , (53) 
and the one body matrix element ((j) + \V PT \(f)~) is proportional to the energy difference 

(0 + |^ PT |0-> ~ ^ + |[<rp,#s-p]|<n ~ - ^ > (54) 

where and e^- are the single-particle energies. In the case of a reflection symmetric (e.g. 
quadrupole deformed) nucleus E + = and E = e^- . The mixing coefficient is therefore 
enhanced only in a few cases when single-particle deformed levels are accidentally close and 
the approximation of Eqs. ( |5B|J5^ ) becomes a crude one [T3|]. In the reflection-asymmetric 
case always [ |!4] , |4Ti| , 



^-^-TO* 1 *^ 1 ' (55) 

where V^d is the reflection-asymmetric part of the single-particle Hamiltonian H s _ p . There- 
fore the admixture coefficient oca is generally enhanced relatively to that in spherical or 
quadrupole deformed nuclei. The splitting within the parity doublet (for K ^ |) can be 



also written as 14,40 



\E + -E\ ~ \E c n p \ , (56) 

exhibiting the reduction of splitting for an odd nucleus in comparison with the splitting of 
parity doublet states of the even-even core. 

Let us look now at the single-particle part of the SM given by Eqs. (^,^,^). The 
"diagonal" single-particle contribution in Eq. ( |49a| ) can be written as 



2a ii (# h |dp|0 i )(a+ i + a. 

V PT \<P-) 



%v+ i) 



« 2(#|O p |^r) pi 1- 0- - <) ■ (57) 

Thus, in general, due to the enhancement in e% the single-particle contribution of the odd 
proton is enhanced relatively to the case of a quadrupole deformed but reflection symmetric 
nucleus. In the proton-odd nuclei the single proton contributes to both terms of the SM 
operator in Eq. (^). These terms are of the same order and partially cancel. In the neutron- 
odd nuclei the single neutron contribution to the SM operator in Eq. (Bj) enters only through 
the corrections to the dipole moment term and therefore is very small. 

As is shown in Ref. (|TJJ the effects of BCS-like pairing are maximal when the Fermi 



level is exactly halfway between the single-particle levels. In this case the simplest model 
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of two opposite parity single-particle levels results in an exactly degenerate parity doublet. 
In realistic calculations the effect of pairing correlations on single-particle levels is revealed 
only in somewhat altered energy splittings between members of parity doublets. In the 
calculations of matrix elements of parity and parity and time reversal violating potentials 
and single-particle parts of electric moments, the effect of pairing is more important. When 
one uses the BCS quasiparticle states instead of single-particle states, the matrix elements 
of the operators V PT , V PY and O p are to be multiplied by pairing factors UfUi + c syra VfVi, 
where c sym depends on the symmetry properties of the operator under the time reversal and 
Hermitian conjugation [2l|,29[. For Hermitian operators, c sym = 1 if an operator is T-odd 



and c sym = —1 if it is T-even p9|| . Thus, for the T-odd V FT , one has the factor UfUi + VfVi, 
whereas for the T- even V PY and a single-particle matrix element of an electric operator 
O p the corresponding factor is UfUi — VfVi. The difference in signs in the expressions for 
pairing correlation factors is related to the fact that matrix element of V PY operator between 
members of a parity doublet goes to zero when the doublet becomes degenerate while the 
matrix element of V PT operator does not vanish p3|,|38|,^4 . 



C. The P- and T- odd interaction 



The P- and T- odd nucleon-nucleon two-body potential can be written in the form pp|f27|] 



as 



W, 



ab 



G 1 

= —^TT- ((Vab<Ta ~ VbaV b )V 5 a (r a - T b ) 

+Vab[ (T a x (T b ]{(p a - pb), 5(r a - r b )} 



(58) 

where G = 10 _5 /m 2 is the Fermi constant (m is the nucleon mass), a,b designate a proton 
or neutron and the curly brackets denote the anticommutation operation. Note, that rj ab are 
in fact effective constants. In this work we use the corresponding effective one-body P- and 
T- odd potential p|,|3,|8 



^|^X>.(v,/ ( ,,)) 



(59) 



where pt(r) = Pof(r) is the nuclear density. We assume |30| , p7| that effective constants r] ab 
and rj are of the same order in magnitude. For the deformed nuclear density we use the 



expansion [^|] 

p(r) = p (r) - fl^^flii - ±-£ff + l(r - RJHofe&VY, 
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(60) 
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where po(r) has the usual Woods-Saxon form, are the nuclear deformation parameters in 
Eq. © andc= T7 ^. 



V. CALCULATIONS AND ESTIMATES 

A. Nuclear structure calculations 

The set of deformation parameters (3i (up to I = 6) was taken from Ref. fl3| where it 
was calculated using a deformed Woods-Saxon single-particle potential and Strutinsky shell 



correction method. From the sets of deformations, calculated in Ref. ||13|| , we have chosen 
those which correspond to the experimental value of K in the g.s. This set is shown in Table 
|. The particle-core wave functions were calculated using the reflection-asymmetric particle 
core model ||14| . We calculated both the intrinsic T-, P- odd moments and admixture coeffi- 



cients using basically the same single-particle potential and shell correction parameters. The 
computer code WSBETA j|2| and the "universal" set of parameters was used in calculations 



involving deformed single-particle Woods-Saxon potential. (A modification of the code was 
made to include explicitly the f3i deformation.) 

To check the stability of the values of admixture coefficients we calculated them also using 
a Nilsson potential with e\ deformations approximately corresponding to the same nuclear 
surface (see e.g. Ref. |[45| ). It is known j|6| that there are some differences in the energies 



and wave functions when calculated with deformed Nilsson or Woods-Saxon potentials in 
the actinide region. Especially the proton levels are different because of the fact that the 
Coulomb term is included in the Woods-Saxon potential but is only simulated in the Nilsson 
potential. Because of that we performed calculation of a - admixtures with Nilsson potential 
only for the neutron-odd nuclei. For the Nilsson potential we used parameters of Ref. |46| 
which are known to produce a good fit to experimental data. See for example the reflection- 
asymmetric calculations of 225 Ra in Ref. ]47| . 



The energy splittings of the core parity doublets - E c and moments of inertia were 
taken from Ref. |34j]. The values of E c we used are given in Table |. The effect of the 



I = 1 deformation (f3\ or ei) in the deformed single-particle potential on the energies and 
admixture coefficients of doublets is small. However this deformation is important for the 
calculation of the shell correction to the intrinsic Schiff moment. 



For the BCS treatment of pairing we used the strength parameter G Pin [48 



1/ N — Z\ 

G p , n = -[g ±gi A ) , (61) 

with g = 19.2 MeV, g x = 7.4 MeV for neutrons and g = 22.0 MeV, g 1 = 8.0 MeV for pro- 
tons. We found that the admixture coefficients, as well as other properties of these nuclei 
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calculated in Ref. |3|], are not sensitive to the pairing strength (or gap parameter A). The 
standard Strutinsky shell correction calculation of the smoothed single-particle occupation 
numbers n was performed using parameters taken from Ref. p7 | . 

The intrinsic SM was calculated using Eqs. ( ^0| , |^T| ) . The "macroscopic" term Sf£ tr was 
computed directly using Eqs. (|5,32) with C = 27 MeV which is the droplet model value 
||32|| . We calculated the shell correction for the first term of the intrinsic SM operator in 
Eq. (||). This correction strongly depends on the proton number and for 86< Z< 91 its 
absolute value decreases with increasing of Z. For Z<92 it has a negative sign relatively 
to Sf^ tr and changes from 33% for 2 f|Rn to 9% for 2 fjAc. In case of proton-odd nuclei we 
calculated also the single proton contribution to the SM. For the g.s. of 223 Fr, 225 Ac and 
229 Pa it amounts to 2-5% of the corresponding values of Sf^ tT . The resulting intrinsic Schiff 
moments of octupole deformed nuclei are in the range (15-28)e fm 3 . We estimate that when 
we allow for reasonable changes in the parameters used, and when other corrections (not 
treated here) are introduced the intrinsic SM will change by 30% at most. Thus we believe, 
that the uncertainties in the evaluated values of the intrinsic SM are of the same order. 

The computed values of intrinsic Schiff moments, admixture coefficients, and resulting 
Schiff moments in the laboratory system as well as the calculated and experimental energy 
splitting for the g.s. parity doublets and calculated parities n p of the intrinsic single-particle 
g.s. are all given in Table 

The main uncertainty in the entire calculation of Schiff moments in the laboratory frame 
arises from the estimate of admixture coefficients a, which are calculated using theoretical 
values of the energy splitting between members of the doublets. In our calculations the 
first state above the Fermi level which has the same value of K as the experimentally 
determined g.s. was chosen to be the g.s. level. In the work of Leander and Chen MM the 
nonadiabatic Coriolis coupling and other refinements were introduced (in some cases also 
adjustments of quasiparticle energies were made), which allowed to describe properties of 
g.s. and excited levels. We did not include couplings between states with different K and 
did not adjust model parameters to fit properties of individual states. In some cases the 
states calculated as the lowest ones above the Fermi surface do not have the experimentally 
determined values of K [[Tj|]. For example the I = | g.s. of 223 Rn which was described in 
34| as arising from Coriolis coupling of K = | and K = | states, is not the lowest state in 



our calculation but is 270 keV above the g.s. Thus it cannot be ruled out that in some cases 
some single-particle states close to the Fermi level and different from the ones we chose, 
fit better the g.s. parity doublet. The absolute values of the admixture coefficients for the 
doublet which was taken as the g.s. are for all nuclei we calculated (except 229 Pa which has 
an exceptionally small energy splitting in the g.s. parity doublet) in the range (1-5) x\0~ 7 r] 
for both reflection asymmetric Woods-Saxon and Nilsson single-particle potentials. If we 
consider also admixture coefficients for the two parity doublets closest to the g.s. parity 
doublet which have the same K the range is (0.15-5) x lO -7 ^. We remark that because the 
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one-body P-, T- odd potential V PT in Eq. (|5^) is proportional to the derivative of the density 
of nucleons and thus is surface-peaked, its matrix elements depend on the behavior of the 
single particle wave functions in the nuclear surface region and hence may vary significantly 
from level to level. 

As one sees in Table some of parities 7r p depart considerably from ±1 meaning that 
these orbitals are strongly parity-mixed, which is one of the reasons for large admixture 
coefficients a. Eq. fl56|) describes well the splittings within the doublets except for the K = ~ 
cases of the 225 Ra and 221 Fr g.s., where there is an additional Coriolis splitting [[14],[38|]. For 
225 Ra, 221 Fr and 225 Ac the experimental energy splitting between members of the g.s. parity 
doublet is well reproduced in our calculation whereas in 223 Fr and 223 Ra nuclei it differs by a 
factor of 2 or 3. For 223 Rn no data on doublets are available. We should remark here that in 
view of Eq. fl5"4|) one can expect some correlation between the matrix element of V PT and the 
energy splitting within the doublet. We expect therefore that the use of experimental energy 
splittings will not necessarily lead to more precise values for the admixture coefficients a. 

The nuclear Schiff moments for the octupole deformed nuclei calculated here are about 
two orders of magnitude larger that those obtained in Refs. P^|J3H|] for isotopes of 129 ' 131 Xe, 
199,201 Hg and 203 > 205 T1. As already mentioned the enhancements we discuss are for nuclei 
with asymmetric shapes and therefore there is also two orders of magnitude enhancement 
with respect to nuclei such as 161 Dy and 237 Np ]18|] which have large quadrupole deformation. 
These nuclei have close to the g.s. levels of the same spin and opposite parity as the g.s., 
but are believed to be reflection symmetric. 

The SMs in even Z nuclei such as Xe and Hg are caused mainly by the polarization of 
protons in the core by the T-, P- odd field of the external nucleon (see Eq. Therefore, 
the SM is proportional to the r] np constant. It was demonstrated in |3(J that this "polariza- 
tion" mechanism gives SM of the same order of magnitude as in nuclei that have a proton 
outside the core. Many-body corrections give also contributions proportional to other T-, 
P- odd constants, e.g. r] nn . In our calculations the interaction of the odd proton or neutron 
with the even-even core and the nuclear SM in the laboratory system are expressed via the 
effective constant r\. 

The experimental data regarding the nuclei we consider are discussed in the reviews 
|,PHHT2]. The 223 Ra and 225 Ra nuclei were also considered in detail in Refs. 0,0. In 



both 221,223 Rn isotopes the I = \ S- s - s P m was determined using laser spectroscopy methods 
5D| , and no data on excited levels are presently available. The computed deformations in 



these two isotopes calculated in Ref. [13| are similar. Using the Coriolis coupling Leander 
and Chen were able to reproduce spectroscopic characteristics of g.s. somewhat better 
for 223 Rn than for 221 Rn. We made a calculation for 223 Rn because its I — | g.s. is easier to 
interpret theoretically Hl3| , |3^| . 

The case of Fr is especially interesting in light of the recent experiments involving trap- 
ping of Fr atoms W% - The latest experimental and theoretical studies of 221 Fr [52] and 223 Fr 
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53| , |54|1 provide strong evidence of intrinsic reflection asymmetry in the g.s. The g.s of 221 Fr 



is K n = \ , I n = | . The assignment of the doublet K n = | + , I n = | + state at 234 keV 
was suggested in ||52|| . The result for the Schiff moment of 221 Fr is smaller than that of 223 Fr 
because of the following factors. Firstly, the factor jr^n m Eq. fl50|) is | for 221 Fr, whereas 
for I = |, K — | g.s. of 223 Fr it is |. Secondly our calculation gives the admixture coefficient 
a for 221 Fr about 3 times less than for 223 Fr. Because of these factors the Schiff moment 
of 221 Fr in the laboratory frame is about 12 times less than that of 223 Fr, although in the 
intrinsic frame the values of Si ntr are roughly equal. 

There is a controversy regarding the g.s. spin and parity doublet in 229 Pa, which is on the 
border of the region of octupole deformed nuclei. Two assignments: a) K = |, I = | g.s. 
and 220 eV energy splitting within the parity doublet [0,0 and b) K = I = \ g.s. and 
unidentified parity partner level [[56T , |1"2"1 were made. In case a) our calculations give using 
the experimental value of the energy splitting the admixture coefficient a = 640 x 10 7 ?7 and 
the Schiff moment in the laboratory frame, S = 230000 x 10~ 8 77 e cm. Note however that in 
Table [TI] the results for 229 Pa are an order of magnitude smaller because a theoretical value 
of 5 keV was used for the energy splitting of the doublet. 



B. Calculation of Atomic Electric Dipole Moments 

Atomic electric dipole moment can be calculated using Eq. (0). However, we do not need 
new complicated numerical calculations to find the EDM of interest. We can use numerous 
calculations for the lighter atomic analogs (Xe, Hg, Cs) and introduce corresponding correc- 
tions taking into account the Z-dependence of the effect to find EDM of heavy atoms (Rn, Ra 
and Fr correspondingly). Indeed, it follows from the atomic calculations that atomic EDM 
in Eq. (|7|) is saturated by the contributions of electrons from the external shells which are 
similar in the analogous atoms (the energy denominators for the transitions from these shells 
are small and radial integrals are large). The expression for the atomic EDM is a product of 
three factors: matrix elements of the radius (fci|r 2 |/c2), energy denominators Ek x — Ek 2 and 
matrix elements of the T-, P- odd nuclear electric potential (ki\ — e(p^\k 2 ). The first two 
factors are determined by the wave functions at large distances and they are the same in anal- 
ogous atoms. This fact is deduced from numerous semi-empirical and computer calculations, 
from experimental data for energy levels and probabilities of electromagnetic transitions, as 
well as from the data on atomic polarizabilities for analogous atoms (the expression for the 
polarizability also contains radial integrals and energy denominators). The matrix elements 
(ki\ — &f^\k2) are determined by the wave function at small distances (more accurately, by 
the gradient of the external electron density at the nucleus) which strongly depends on the 
nuclear charge. However, this dependence was calculated analytically [p~S|| . The contribution 
of the Schiff moment is determined by the matrix element between sy2 and py2 (o r S1/2 and 
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P3/2) electron orbitals which is proportional to SZ 2 R 1 / 2 (or SZ 2 R 3 / 2 ). The number of p 3 / 2 
states is two times larger then the number of P1/2, therefore we will need linear combination 
of the relativistic factors R sp = (R1/2 + 2i?3/2)/3. 
The relativistic factors R1/2 and R3/2 are given by 

4 7l/2 (2ZR \ 2 ^/*- 2 



R3/2 



[T(2 7 i/2 + I)] 2 V a B 

48 (2ZR \ ti/2+73/2-3 

r(2 7 i /2 + i)r(2 73 /2 + 1) 



(62) 



where 7j - = [(J + |) 2 — (Za) 2 ] 1 / 2 , a s is the Bohr radius and i? is the nuclear radius. In light 
atoms R sp ~ 1, in heavy atoms R sp ~ 10. Thus, we have simple estimates of the electric 
dipole moments of Ra, Rn and Fr induced by the nuclear Schiff moments: 



dat(Ra) = dat(Hg) 
4t(Rn) = 4t(Xe) 



{SZ 2 R 

(SZ 2 R sp ) Kg 

(S Z 2 R sp ) Rn 
(SZ 2 R, 



spJXe 



4t(Fr) = d a t(Cs) f q f 2 l sp \^ , (63) 

where only the nuclear Schiff moment contributions to atomic EDMs rf at of Hg, Xe and Cs 
are taken into account. The atomic structure ratios here are 
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(Z 2 R sp )x e 



7.7 . (64) 



The calculation of EDM of Hg, Xe and Cs have been done in Refs. |lqj57| , |30[] . The results 
of our calculations for Ra, Rn and Fr are presented in the Table |I|. 



VI. SUMMARY 

In this paper we studied of the T-odd, P- odd electric moments in heavy nuclei with 
intrinsic reflection asymmetry and induced electric dipole moments in corresponding atoms. 
We presented a detailed theory of the collective T-, P- odd electric moments in reflection 
asymmetric odd- mass nuclei, in particular the Schiff moment. We employed the two-fluid 
liquid-drop model, particle plus core model, and used the results of Nilsson-Strutinsky mean 
field calculations for intrinsic reflection asymmetric nuclear shapes. Various corrections for 
nuclear T-odd, P-odd electric moments were evaluated. 
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In the calculations of induced atomic electric dipole moments we employed the scaling 
relations between such moments in heavy atoms and their lighter analogs and used the 
results of the calculations for the latter to find the corresponding moments in heavy atoms. 
We studied the cases of all the heavy reflection-asymmetric odd-mass nuclei for which there 
is evidence of intrinsic octupole deformation in the ground state and which are relatively 
long-lived, so their atoms could be suitable for experiment. 

The results can be summarized in the following form: 

a) In a reflection asymmetric nucleus which has odd mass number or is odd-odd, enhanced 
collective T-odd, P-odd electric moments appear, if T-odd, P-odd terms are present in the 
nuclear Hamiltonian. 

b) The T-odd, P-odd Schiff moments in heavy nuclei with intrinsic reflection asymmetry 
are typically enhanced by more than two orders of magnitude in comparison with reflection- 
symmetric deformed or spherical nuclei. 

c) Due to the atomic structure effects, atomic electric dipole moments in heavy atoms are 
enhanced, compared to the lighter analogs. For atoms of nuclei with Z around 90, the 
atomic enhancement is of about 8 times, in comparison with analog atoms with Z around 
55. This enhancement factor is about 2 compared to analogs with Z around 80. 

d) The atomic electric dipole moments, induced by T-odd, P-odd hadron-hadron interaction 
in the nuclei studied are typically enhanced 400-1000 times in comparison with Hg and Xe 
nuclei, for which the best experimental upper limits on atomic electric dipole moments are 
obtained. These findings may open new experimental possibilities of studying time reversal 
violation. 
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APPENDIX A: SCREENED T-, P- ODD ELECTROSTATIC POTENTIAL OF A 
NUCLEUS AND THE SCHIFF THEOREM 

The Hamiltonian of an atom placed in a homogeneous external electric field E is 
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H = J2(K l - etp (Ri) - eR t E ) 



+ E 



i>k 




-dE t 
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^)-/^. (Al) 
J \Ki — r\ 

Here K, t and Ri are the kinetic energy and coordinate of the electron, (p Q (Ri) is the electro- 
static nuclear potential and d is the nuclear dipole moment. Let us add to H an auxiliary 
term 

V = dE -^-J2 dV m {Ri) • (A2) 
ez i 

It is easy to demonstrate that in the linear approximation in d the interaction V does not 
produce any energy shift, (V) = 0. Indeed 



i 

m 



[]T Pl , H] = -e Vi<pa{Ri) + ZeE . (A3) 



We have taken into account that the total electron momentum J2iPi commutes with the 
electron-electron interaction term. Using Eq. ( |A2| ) and (n\[H, J2iPi]\ n ) ~ (En — E n ) = we 
obtain (V) = dE - -^eZdE = 0. 

To find an electric dipole moment one needs to measure a linear energy shift in an external 
electric field. Since V does not contribute to this shift we can add it to the Hamiltonian 



H = H + V = Y^(Ki - ecfiRi) - eR t E Q 

i 

+ E 



e 2 



i>k I-*** _ R k\ 

<p(Ri) = <f (Ri) + —dVtfo{Ri) ■ (A4) 
eZ 

Note, that the Hamiltonian H does not contain the direct interaction dEo between the 
nuclear electric dipole moment and external field (Schiff theorem [32,231). The dipole term 
is also canceled out in the multipole expansion of <p(Ri). 
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TABLES 

TABLE I. Intrinsic g.s. deformations and energy splittings between opposite parity core states. 





223 R a 


225d„ 
IX cl 


223p n 
XVI 1 


221 Fr 


223 -p 


225 A P 


229p„ 


Pi 


n 1 9^ 


n i a ^ 

U. L^tO 


n 1 9Q 
u. izy 


U. 1UO 


199 


U. loo 


U.l I D 


& 


0.100 


0.099 


0.081 


0.100 


0.090 


0.104 


0.082 


Pa 


0.076 


0.082 


0.078 


0.069 


0.076 


0.078 


0.093 


& 


0.042 


0.035 


0.024 


0.045 


0.033 


0.038 


0.020 


As 


0.018 


0.016 


0.023 


0.020 


0.022 


0.013 


0.015 


E c (keV) 


212 


221 


213 


305 


212 


206 


333 



TABLE II. Admixture coefficients a (absolute values) and theoretical energy splitting be- 
tween the g.s. doublet levels AE = E~ — E + , parities of the intrinsic (reflection-asymmetric) sin- 
gle-particle g.s. calculated using the Woods-Saxon ( W) and Nilsson (Nl) potentials, experimental 
energy splitting, intrinsic Schiff moments and Schiff moments (with Woods-Saxon potential) as well 
as induced atomic dipole moments. The values for 199 Hg, 129 Xe and 133 Cs from Refs. ]30| , |5"7 ,18| 
are given for comparison. 
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225 Ra 


223 Rn 


221 Fr 


223JY 
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a(WS) (10 7 rj) 
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A£(WS) (keV) 
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75. 
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5. 






7Tp(WS) 


0.81 


-0.02 


0.17 


-0.55 


-0.34 


-0.35 


0.01 






a(Nl) (10 7 77) 


2. 


5. 


2. 














AE(N1) (keV) 


171. 


55. 


137. 














AE cxp (keV) 


50.2 


55.2 




234. 


160.5 


40.1 


0.22 






S'intr (e fm 3 ) 


24 


24 


15 


21 


20 


28 


25 






S (10 8 7] e fm 3 ) 


400 


300 


1000 


43 


500 


900 


1.2xl0 4 -1.4 


1.75 


3 


d(at) (10 25 7? e cm) 


2700 


2100 
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240 


2800 




5.6 


0.47 


2.2 
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